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The full-potential linearized augmented plane wave (FP-LAPW) method with the generalized 

gradient approximation (GGA) for the exchange-correlation potential has been applied to predict 

the elastic properties of quaternary oxyarsenide LaOFeAs - the basic phase for the newly discovered 

00 , 26-52K superconductors. The optimized lattice parameters, independent elastic constants (Cij), 

f^ ■ bulk modulus B, compressibility (3, and shear modulus G are evaluated and discussed. For the 

("V^ ' first time the numerical estimates of the elastic parameters for poly crystalline LaOFeAs ceramics 

Cn \ are performed. Our analysis shows that LaOFeAs belong to the anisotropic mechanically stable 

}_( . systems; a parameter limiting the mechanical stability of this material is the shear modulus. In 

Qh' addition, LaOFeAs is a soft material (B ~ 98 GPa) with high compressibility {13 ~ 0.0102 1/GPa); 

■^\ I lays at the border of a brittle behavior and adopts a high ionic contribution in intra-atomic bonding. 
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The elastic properties are of great importance for the 
material science of superconductors; on the other hand, 
the elastic constants can be linked to such important 
physical parameters of superconductors, as the Debyc 
temperature Qo and the electron-phonon coupling 
constant A [l|. In the past years, the elastic proper- 
ties of various superconducting materials have been 
studied extensively, and some correlations between the 
superconducting critical temperature Tc and mechanical 
parameters have been discussed. So, Hirsh [2] has 
suggested that high TcS are associated with low values 
of the bulk modules B, i.e. high compressibility (3. 
Really, for the many groups of superconductors with 
enhanced critical temperatures (such as YBCO, MgB2, 
MgCNia, boride-carbides and carbide- halides of the rare 
earth metals, REM2B2C and RE2C2X2 etc) their bulk 
modules does not exceed B < 200 GPa (/3 > 0.005 
1/GPa) [H, i, i, i, i, i, Q- On the other hand, the 
superconducting transition (up to Tc ^ IIK) was found 
for such hard and incompressible material as boron - 
doped diamond [a, Q- 

Except scientific interest, mechanical properties appear 
extremely important for technology and various ad- 
vanced applications of superconducting materials 10']. 
In February 2008 the new layered superconductor: 
electron doped quaternary oxyarsenide LaOFeAs 
(LaOi_:rF2,FeAs,_^x~0.05-0.12) with T^ at about 26K 
was reported [111 ], and this surprising discovery has 
stimulated much activity for search of the new supercon- 
ducting materials based on this phase. It was found that 
also hole-doped oxyarsenide LaOFeAs (Lai-xSr^jFeAs, 
by partial substitution of La by Sr) adopts supercon- 
ductivity with Tc at about 25K 12]. Moreover by 
replacing of La by others rare earth elements it seems to 
be possible to achieve the further enhancing of Tc (41K 
- for CeOi_^F^FeAs [H, 43K - for SmOi-^^F^rFeAs 
lli_5p-52K - for PrOi_^F^FeAs and NdOi-^^F^rFeAs 

Hill). 

Alongside with the intense efforts in synthesis of new 
doped quaternary oxypnictides and examining of su- 



perconductivity parameters, some of their physical 
properties have been investigated. So, today the inter- 
play between superconductivity and spin fluctuations, 
band structure, Fermi surface, and phonon spectra 
for some quaternary oxypnictides are discussed at the 
theoretical level [11 El, d, S [SI, [H, IH, [H; the 
photoemission spectra of SmOo.ssFo.isFeAs are ob- 
tained [2^] and pressure exp eriments for La0i_a;F2,FeAs 
(x=0.11) are performed 129]. 

In this Report we focus on the elastic properties of 
monocrystalline quaternary oxyarsenide LaOFeAs as 
the basic phase for the newly discovered 26-52K su- 
perconductors. In addition the numerical estimates of 
the elastic parameters of the polycrystalline LaOFeAs 
ceramics, which may be important for the future appli- 
cations of new superconducting materials are performed 
for the first time. 

The oxyarsenide LaOFeAs has a tetragonaj_(ZrCuSiAs- 
type structure, space group P4/nmm[27|) layered 
structure, where Fe atoms are arrayed on a square 
lattice in the form of edge-shared tetrahedra FeAs4. 
Each (Fe-As)^ layer is sandwiched between (La-0)+ 
layers. 

Our calculations based on density functional theory 
(DFT) with generalized gradient approximation (GGA) 
in PBE form [23] for the exchange-correlation potential 
and were performed by means of the full potential 
method with mixed basis APW-l-lo (LAPW) imple- 
mented in the WIEN2k suite of programs [29|. The 
calculations are performed with full lattice optimizations 
(lattice parameters and atomic coordinates) ; the self- 
consistent calculations were considered to be converged 
when the difference in the total energy of the crystal 
did not exceed 0.001 mRy as calculated at consecutive 
steps, and the nonmagnetic state was treated. Other 
details are described in Ref [6|]. 

Firstly, the equilibrium lattice constants (a and c) were 
evaluated - in reasonable agreement with available data, 
Table 1. 
Secondly, the values of six independent elastic constants 



TABLE I: Calculated lattice constants (a, c, in A), cell vol- 
ume (Vo, in A^) elastic constants {Cij, in GPa), bulk modulus 
(B, in GPa), shear modulus (G, in GPa) and B/G ratio for 
tetragonal monocrystalline LaOFeAs as well as some elastic 
parameters for polycrystalline LaOFeAs ceramics as obtained 
in the Voigt-Reuss-Hill approximation: bulk modulus (Bvrh, 
in GPa), compressibility (/3 vrh, in GPa~^), shear modulus 
{GvRH, in GPa), Youngs modulus {Yvrh, in GPa) and Pois- 
sons ratio {v). 



a 


4.033 (4.036 [11]; 4.037 [30])" 


Bv(Bfl)'' 


98.5 (97.2) 


c 


8.684 (8.739 [11]; 8.629 [30])" 


GviGR)" 


56.5 (55.9) 


Vo 


141.2 


Hv /Gv 


1.74 


Gil 


191.9 


BvflH 


97.9 


Gl2 


55.9 


/3 VRH 


0.01022 


Cl3 


61.6 


GvRH 


56.2 


G33 


144.8 


Yvrh 


141.5 


G44 


44.1 


V 


0.259 


Gee 


77.9 







" in parentheses the available experimental (Ref . [Till ) and calcu- 
lated (Ref . ^31 ) data are given 
'' in Voigt (and Reuss) approximations 



(Qj-; namely Cn, C12, C13, C33, C44 and Cee) for 
LaOFeAs were estimated by calculating the stress 
tensors on applying different deformations given on the 
equilibrium lattice of tetragonal unit cell and determin- 
ing the dependence of the resulting change in energy 
on the deformation. Table L All these elastic constants 
are positive and satisfy the well-known Born's criteria 
for mechanically stable tetragonal crystal: Cn > 0, 
C33 > 0, C44 > 0, C66 > 0, (Cn - C12) > 0, (Cn + 
C33 - 2C13) > and 2(Cn + C12) + C33 + 4Ci3 > 0. 
The calculated anisotropic factor A= 2C44/(Cii-Ci2) is 
0.65, i.e. LaOFeAs phase is mechanically anisotropic, 
because a value of A = 1 represents completely elastic 
isotropy, while values smaller or greater than 1 measure 
the degree of elastic anisotropy. 

The calculated elastic constants allow us to obtain the 
macroscopic mechanical parameters of LaOFeAs, namely 
bulk modulus (B) and shear modulus (G) - in two main 
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approximations: Voigt (V) [3]| and Reuss (R) [32], in 
the following forms: 



By = -{2(Cn + C12) + C33 + 4C13} 



Gv = —{M + 3Cn - 3Ci2 + I2C44 + eCge} 



M 



Gr = 15{ 



18Bv 



' +^}-^ 



C^ Cn — C12 C44 Cee 



where C^ = (Cn + 012)033 - 2C13 and M = On + C12 
-I- 2C33 - 40 13. The results obtained are summarized in 
Table L 

Next, as LaOFeAs samples are usually prepared and 
investigated as polycrystalline ceramics [11,12,25,26] in 
the form of aggregated mixtures of microcrystallites 
with a random orientation, it is useful to estimate 
the corresponding parameters for the polycrystalline 
material from the elastic constants of the single crystal. 
To this aim we utilize the Voigt-Reuss-Hill (VRH) 
approximation. In this approach, according to Hill[33], 
the Voigt and Reuss averages are limits and the actual 
effective modules for polycrystals could be approximated 
by the arithmetic mean of these two bounds. Then, one 
can calculate the averaged compressibility (/3 vrh = 
l/By/jff), Young modulus (Yvrh) and from the Bvrh, 
GvRH and Yvrh one can evaluate the Poissons ratio 
(z^). All these parameters are also listed in Table L 
Certainly all of these estimations are performed in limit 
of zero porosity of LaOFeAs ceramics. 
From our results we see that for LaOFeAs Bvrh > 
GvRH', this implies that a parameter limiting the me- 
chanical stability of this material is the shear modulus 

GvRH ■ 

The bulk modulus for LaOFcAs is rather small (at about 
98 GPa), and is less than, for example, the bulk modules 
for other known superconducting species: MgB2 (^ 
122-161 GPa [34]), MgCNi3 (-172-210 GPa [35]), YBCO 
(- 200 GPa), and YNi2B2C (-200 GPa [3]). Thus, 
compared with other superconductors, LaOFeAs is the 
soft material. 

According to criterion [36] , the material is brittle if the 
B/G ratio is less than 1.75. In our case, for LaOFeAs 
B/G — 1.74, this suggests that this material lays at the 
border of a brittle behavior. 

Finally, the values of the Poisson ratio (v) for covalent 
materials are small {v = 0.1), whereas for ionic materials 
a typical value of v is 0.25 [37]. In our case the value 
of h' for LaOFeAs is at about 0.26, i.e. a higher ionic 
contribution in intra-atomic bonding for LaOFeAs phase 
should be assumed. Besides, for covalent and ionic 
materials, the typical relations between bulk and shear 
modulus are G ~ I.IB and G~ 0.6B, respectively. 
For our case the calculated value of Gy/By is 0.58, 
also indicating that the ionic bonding is suitable for 
LaOFeAs. 

In summary, we have performed FLAPW-GGA calcu- 
lations to obtain the elastic properties for mono- and 
polycrystalline oxyarsenide LaOFeAs as a basic phase 
for the newly discovered 26-52K superconductors. Our 
analysis shows that LaOFeAs belong to the anisotropic 
mechanically stable systems; a parameter limiting 
the mechanical stability of this material is the shear 
modulus. In addition, LaOFeAs is a soft material 
(B-98 GPa) with high compressibility {f3 - 0.0102 
1/GPa); lays at the border of a brittle behavior and 
adopts a high ionic contribution in intra-atomic bonding. 
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